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■^-j- ■ Abstract 

A quantum mechanical system describing bosons in one space dimension with a kinetic 
■ energy of arbitrary order in derivatives and a delta function interaction is studied. Exact 

wavefunctions for an arbitrary number of particles and order of derivative are constructed. 
| Also, equations determining the spectrum of eigenvalues are found. 

PACS numbers: 02.30Jr, 05.30Jp 

Physical systems whose equations contains derivatives higher than two arise in a number of 
contexts. In classical mechanics, for example, the Korteweq-deVries equation, which describes 
shallow water waves, contains a third order derivative term [HI. Higher derivative theories of 
classical gravity, first introduced by Weyl |2j , have some attractive cosmological properties, such 
as gravity driven expansion ||. A model of quantum gravity with higher derivative terms has 
been shown to be renormalizable H and asympotically free In particle physics, the Nambu- 
Jona-Lasinio ||, Schwinger [0, and Skyrme || models with higher derivative terms added have 
been considered. Also, in the superfield formulation of supersymmetric field theories, higher 
derivative terms arise naturally 0. 

In this work, a family of higher derivative generalizations of Schrodinger's equation is studied. 
The system describes bosons where the "kinetic energy" is a polynomial of order n in derivatives 
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and the interaction is a delta function potential. A particular choice of polynomial corresponds 
to the much studied delta function bose gas (or quantum nonlinear Schrodinger equation) [[II], Q |. 
The main result here is that for an arbitrary polynomial in derivatives and arbitrary number of 
particles, exact eigenfunctions for the system may be constructed. 

Formulating the system under study, consider the eigenvalue problem 

H^ = E^, (1) 

where the hamiltonian H is the partial differential operator 

N n 

H = EE^K-) 1 + E 2c <K^ - x k), (2) 
j=i 1=1 j<k 

and ffr(xi, ...,xn) is a symmetric function of the Xj. N is the number of bosons, n is an arbi- 
trary natural number, the ai are arbitrary real coefficients, and c is a real interaction constant. 
Equation (j^) is a generalization of Schrodinger's equation describing bosons interacting via a 
delta function where the usual 2nd order kinetic energy has been replaced with a polynomial 
containing higher order derivatives. 

The hamiltonian (fj) is translationally invariant, thus the momentum operator 

N 

P = -iJ2d Xj (3) 
i=i 

commutes with the hamiltonian (0). 

If a 2 = 1, and all other ai = 0, then (H) is called the 5-function bose gas, and the eigenfunctions 
for arbitrary iV were first found by Lieb and Liniger [ID|. The method of solution has come 



to be known as coordinate Bethe ansatz, after Hans Bethe, who first used the idea to find the 
eigenstates of the Heisenberg spin chain [12], and has since been applied to many 1 + 1 dimensional 
quantum and 2 dimensional statistical mechanical systems fL3| . 

Proceeding now with constructing the the eigenfunctions of (0), first consider a simpler model 
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where a n = 1 and all other ai are equal to zero 

N 

H = EH<^r + E 2cS ( x 3 - x k) =T n + V. (4) 

3=1 j<k 

The formulae for the eigenf unctions of this model are simpler and it will be seen that the solution 
for the general model (0) can be easily constructed knowing the solution of (f|). 

Next, note that for coordinate values where no two x^s coincide, Xj ^ Xf., the interaction 
term V is zero and the eigenvalue equation (JIJ) reduces to 

N 

1 E(-id x .) n * = E*. (5) 

3=1 

A "basic" solution of this equation is 
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e(x 1 ,..,x N ;k 1 ,...,k N ) = expi(E^j), (6) 

j=i 

where the kj are a set of quasimomentum parameters. The energy and momentum of this state 
are 

N N 

£ = P = Y,ky (7) 

3=1 3=1 

The effects of the interaction term V may be rewritten as a set of conditions which $ must 
satisfy when crossing a boundary Xj — Xf.. To obtain these conditions, integrate equation (|T|) 
across a boundary Xj = x k . 

lim / d(x k - Xj)H^ = lim f d(x k - Xj)E^ (8) 

The RHS is zero, leaving only the kinetic and interaction terms. Calculating the integral 



lim/ d(x k - Xj )(T n + V)V, (9) 

yields the equation 

' T E (f) (! + (-l) n "0((^ + d X] )\d Xk - <),,)>■ ' M>).,, ,, + c*(z* = x,) = 0, (10) 



yn- 

1=0 



where (jj is the usual binomial coefficient. Thus, the orignal problem of solving (|I]) has been 
recast as that of solving the noninteracting equation (^) subject to the conditions (|rU|). 

To simplify the task of solving (||) and (|i0|) , note that since this system describes bosons on 
a line, it is sufficient to construct the eigenfunction only in an ordered sector, 

XQl < %Q2 < ••• < XQn-i < X Qni ( n ) 

where Q is a permutation of N elements. The value of the wavefunction for other coordinate 
orderings can be obtained by symmetric continuation. For convenience, we set Q — 1. 
Next, consider the following superposition of the solutions (|6]) 

N 

Yl <t>(ki,...,k N ;P)expi(J2k P .x j ), (12) 

P£S N j=l 

where ki,...,k^ are quasimomenta, P is a permutation of N elements, and <f)(kx, k N ; P) are 
N\ undetermined phase functions of the quasimomenta and P. The energy and momentum of 
this state are given by (^). 

Now it will be determined how to choose the <p{P) so that the boundary condition (|10D is 
satisfied across an intersection Xj = Xj+%. Consider the partial sum 

N N 

§ = 0(P)expi(E^i) +<j>{P')e W i{Y,kp'X j ), (13) 

3=1 3=1 

where P' differs from P by the exchange of kp. and kp j+l . If this partial sum satisfies the 
boundary conditions (|T0D, then so does the the entire sum ( p!2|) since given a pair j,j + 1, then 
TV! permutations may be split into ^ pairs of the type in expression (|T3"|). Now requiring that 
the partial sum (|T3|) satisfies (|T0|) , one obtains a set of conditions on the phases 4>{P) 

E (") (1 + (-ir'X^ + k Pj+1 ) l (k Pj+1 - k P] Y- l ~\<P{P) - 0(P')) + c(0(P) + 0(P')) = o. 

(14) 
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To write this in a more transparent form, define the polynomial 

p n (k, *0 = air E (") ( x + WW + - k ) n ' l ~ x - ( 15 ) 

Note that this polynomial is of degree n — 1 and odd under the interchange and fc'. Rewrit- 
ing fll4]) in terms of the ratio of phases, one obtains 

<f,(P') _ Pn(k Pi ,k Pi „)+ic _ 

m ~ p,(kp„k P , tt) ~ t c = e " (kp " kp ^ (16) 

These equations have the interpretation as the phase accrued when two particles scatter off one 
another. Note that the scattering phase has the property l/9 n (k, k') = n (k', k). 
Lemma:A solution to the conditions ( |TED is 

<!>{P)= IT Uk Pl ,k Pm ). (17) 

/>m,P ; <Pfc 

Thus, the result for the eigenfunctions of (f|) is 

N 

*= E <P(ki,-,k N ;P)expiC£kp j x Qj ), (18) 
PeSV i=i 

with 0(P) chosen according to flT7|), locally satisfies the eigenvalue equation (JI|) for all natural n 
and N. The eigenvalues of if and P are given by (|7|). 
Properties of the solution: 

1. The wavefunction \1> is continuous for all Xj, and the first derivative d Xj ^(x) is discontinuous 
when any two x's coincide 

2. ty(kj = ki) = for all j, k, thus a Pauli-type exclusion principle holds for the quasimomenta 
parameters. 

3. Explicit formulae for the first few 6 n are 

n = l: 0i = -1 (19) 



Uk- - if 2 ) + ic 
« = 3: "3= |^_^_ ic (21) 

_ \{7k 3 + 3A; 2 *;' - 3£;A;' 2 - 7k' 3 ) + ic 
n = 4: 0A = \{7k 3 + 3k 2 k' - 3kk' 2 - 7k' 3 ) - ic ^ 

|(15fc 4 + 10k 3 k' - 10kk' 3 - 15A;' 4 ) + ic 
' 5 |(15A; 4 + 10fc 3 A;' - 10/fc/fc' 3 - 15/fc' 4 ) - ic 1 ' 

The case n = 1 corresponds to free fermions with a linear energy spectrum. 

Next, to construct the solution of the general model @, with arbitrary a/, note that since the 
differential equation (P is linear, all polynomials corresponding to different values of I simply 
add. So, stating the main result: 

Theorem: the eigenfunctions of (fj) are of the form 

N 

*= Yl &(ki,...,k N ;P)expiQ2kp j x Qj ), (24) 

where 



HP)= II ®n(kp v k P J, (25) 

l>m,Pi<P h 

and the two body scattering phase is 

n 

Y,aiPi(k,k') + ic 
G n (k, k') = ^ . (26) 

^2aiPi(k, k') - ic 
i=i 

n obeys the properties 1/Q(k,k') = Q(k',k) and Q(k,k) = — 1. The latter implies that the 
Pauli exclusion principle for spectral parameters holds for the general model. The eigenvalues of 
H and P are given by 

N n N 
3=1 1=1 j=l 
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Thus far, the boundary conditions on the partial differential equation ([!]) acts have been left 
unspecified. The next step is to impose that the wavefunction \1/ is periodic in a box of length 
L, ie., require that 

• /-) *!•'-,) ./ 1 V. (28) 

This condition imposes the constraint that the parameters fcj, k N must satisfy a set of coupled 
equations 

N 

e ¥ = -Il8»fey j = l,...,N. (29) 
i=i 

An interpretation of these equations is as the phase acquired by particle j as it traverses around 
the box back to its starting point, scattering off the other particles along the way. 

To understand the ground state, excitation spectrum, and other physical characteristics of 
the hamiltonian (H), one needs to analyze the solutions of the equations The properties of 
equations ( ^9l) can differ essentially from those of equations found for other 1 + 1 dimensional 
quantum systems (e.g. those of fl3fl ). In particular, for complex solutions, the imaginary part of 
a solution can have a nontrivial functional dependence on the real part. The properties of these 
complex solutions need to be understood in order to study physical properties which depend on 
the spectrum of eigenvalues. 

The author thanks Minoru Takahashi for useful discussions and the Japanese Society for the 
Promotion of Science for financial support. 
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